//mumax 3.10 [windows_amd64 go1.11.6(gc) CUDA-11.0]
//GPU info: NVIDIA GeForce RTX 4080(16375MB), CUDA Driver 12.0, cc=8.9, using cc=80 PTX
//(c) Arne Vansteenkiste, Dynamat LAB, Ghent University, Belgium
//This is free software without any warranty. See license.txt
//********************************************************************//
//  If you use mumax in any work or publication,                      //
//  we kindly ask you to cite the references in references.bib        //
//********************************************************************//
//output directory: .\Scripts\Current\SAF_Sky_Breathing_Current_RKKY_-0.00030_50mT.out/
cellsizeXY := 1.5625e-9
cellsizeZ := 0.4e-9
Nx := 256
Ny := Nx
SetPBC(5, 5, 0)
SetGridSize(Nx, Ny, 4)
SetCellSize(cellsizeXY, cellsizeXY, cellsizeZ)
Ms1i := 5.80e+05
Ms2i := 580e3
Msat_value := (Ms1i + Ms2i) / 2
Ku1value := 0.7e6
Msat = Msat_value
Aex = 15e-12
Dind = 3e-3
Ku1 = Ku1value
alpha = 0.01
AnisU = vector(0, 0, 1)
B_ext = vector(0, 0, 0)
MagLayers := Layer(0)
for i := 0; i < 2; i++ {
	MagLayers = MagLayers.add(Layer(2 * i))
	DefRegion(i, Layer(2*i))
}
SetGeom(MagLayers)
Msat.SetRegion(0, Ms1i)
Msat.SetRegion(1, Ms2i)
RKKYAFM := -0.00030
B_ext = vector(0, 0, 0)
TableAdd(B_ext)
TableAdd(CropLayer(m, 0))
TableAdd(CropLayer(m, 1))
TableAdd(CropLayer(m, 2))
TableAdd(CropLayer(m, 3))
TableAdd(m.region(0))
TableAdd(m.region(1))
TableAddVar(RKKYAFM, "RKKY", "j/m^2")
B_AFM1 := Masked(Mul(Const(RKKYAFM/(cellsizeZ*Ms1i)), shifted(m, 0, 0, -2)), layer(0))
B_AFM2 := Masked(Mul(Const(RKKYAFM/(cellsizeZ*Ms2i)), shifted(m, 0, 0, 2)), layer(2))
B_AFM := Add(B_AFM1, B_AFM2)
Edens_AFM := Mul(Msat, Dot(Mul(Const(-1), m), B_AFM))
AddFieldTerm(B_AFM)
AddEdensTerm(Edens_AFM)
for i := 0; i < 2; i++ {
	m.SetRegion(i, NeelSkyrmion(pow(-1, i), pow(-1, i+1)).scale(3, 3, 1))
}
saveas(m, "initialised_m")
B_ext = vector(0, 0, 50e-3)
relax()
saveas(m, "relaxed_m")
lambda = 1
Pol = 0.56
epsilonprime = 0
fixedlayer = vector(0, 0, 1)
jc := 1e9
f := 100e9
J = vector(0, 0, jc*sinc(2*pi*f*t))
tableautosave(2e-12)
Autosave(m, .5e-9)
run(5e-9)
//********************************************************************//
//Please cite the following references, relevant for your simulation. //
//See bibtex file in output folder for justification.                 //
//********************************************************************//
//   * Mulkers et al., Phys. Rev. B 95, 144401 (2017).
//   * Vansteenkiste et al., AIP Adv. 4, 107133 (2014).
